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Abstract 

We study how fluctuations of the black hole geometry affect the properties of 
Hawking radiation. Even though we treat the fluctuations classically, we believe 
that the results so obtained indicate what might be the effects induced by quantum 
fluctuations in a self consistent treatment. To characterize the fluctuations, we use 
the model introduced by York in which they are described by an advanced Vaidya 
metric with a fluctuating mass. Under the assumption of spherical symmetry, we 
solve the equation of null outgoing rays. Then, by neglecting the greybody factor, we 
calculate the late time corrections to the s-wave contributions of the energy flux and 
the asymptotic spectrum. We find three kind of modifications. Firstly, the energy 
flux fluctuates around its average value with amplitudes and frequencies determined 
by those of the metric fluctuations. Secondly, this average value receives two positive 
contributions one of which can be reinterpreted as due to the 'renormalisation' of 
the surface gravity induced by the metric fluctuations. Finally, the asymptotic 
spectrum is modified by the addition of terms containing thermal factors in which 
the frequency of the metric fluctuations acts as a chemical potential. 
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1 Introduction 



25 years have passed since Hawking's theoretical discovery of quantum radiation by black 
holes Since then, many aspects of this phenomenon have been investigated. First, the 
mean value of the energy-momentum tensor has received much attention since it provides 
the source of the so-called semi-classical Einstein equations (see e.g. 0, |3| and references 
therein). Hopefully, the solutions of these equations should govern the mean evolution of 
the evaporating geometry. More recently, more quantum mechanical questions which raise 
doubts concerning the validity of this semi-classical evolution have also received attention. 
In particular, the controversial role of arbitrarily large ('transplanckian') frequencies of 



vacuum fluctuations B- [10] and the gravitational back reaction due to a specific quantum 



TT| , |12fl have been discussed. 

In this paper, we shall consider another aspect: We study how the fluctuations of the 
black hole horizon geometry might affect the properties of Hawking radiation. To describe 
these fluctuations quantum mechanically and to determine their effects on Hawking radi- 
ation requires full quantum gravity. Besides the 'spontaneous' metric fluctuations there 
also exist so-called 'induced' metric fluctuations, which are generated by quantum fluc- 
tuations of all other fields interacting with the gravitational one. In the regime when the 
'induced' metric fluctuations are dominating, a consistent way to describe black hole fluc- 
tuations and back reaction is to use the stochastic semiclassical theory of gravity based on 



the Schwinger-Keldysh effective action [O, [l|J and the Feynman- Vernon influence func- 



tional [15, [L6| methods. In stochastic gravity the semiclassical Einstein equations are 
generalized to Einstein-Langevin equations which contain stochastic stress-energy tensor 



describing metric fluctuations induced by quantized fields. (For recent review see [|T? 



20H and references therein.) 

The study of the effects connected with black hole fluctuations is a technically very 
complicated problem. Only some preliminary work has been done in this direction till 
now. Under these conditions it is natural to study simplified models. In particular, it 
is not unreasonable to hope that the main properties of the Hawking radiation modified 
by metric fluctuations can be extracted from a much simpler framework in which the 
fluctuations of the metric are treated classically. 

The model we shall use is inspired by that proposed by York ||21|| . In that model, the 
fluctuating geometry near the horizon of the black hole is represented by a Vaidya-type 
metric with a fluctuating mass. The spectrum of these fluctuations is characterized by the 
zero point fluctuations of quantum fields. In this paper, we further simplify this model 
by considering only spherically symmetric fluctuations and by neglecting the scattering 
by the gravitational potential which occurs in the 4-dimensional Dalembertian. Then we 
determine how these fluctuations modify the energy flux and the asymptotic spectrum of 
s- waves. 

The paper is organized as follows. York's model is described in Section 2. Section 
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3 contains a perturbation analysis of the equation of radial null ray propagation in the 
fluctuating geometry. The solution of this equation is obtained in Section 4 and used to 
obtain the modified energy flux in Section 5 and the spectrum in Section 6. The results 
are discussed in Section 7. In our work we use dimensionless units where G — c = h — 1 



and the sign conventions of [22 



2 Model 

In order to study the influence of metric fluctuations on Hawking radiation, we consider 
a simplified version of the model proposed by York [ETJ. In his model, the fluctuations of 



the black hole geometry are approximated by an incoming Vaidya metric with a fluctu- 
ating mass. The fluctuating part of the mass function can be decomposed into spherical 
harmonics. Upon quantizing the gravitational field, only components with i > 2 are im- 
portant. However quantum fluctuations of matter fields may induce fluctuations of the 
geometry with all i. 

In what follows we shall only consider spherical modes of fluctuations. Therefore, the 
metric for a neutral non-rotating black hole can be taken of the form 

ds 2 = -A dv 2 + 2dv dr + r 2 dS\ , (2.1) 

where dS 2 is the metric of a unit 2-sphere and 

2m 

A = l , (2.2) 

r 

m = m(v) = M[l + n(v)]-&(v), (2.3) 



n{v) = /Xosin(u;t>) . (2.4) 

This is the standard Vaidya metric in advanced time coordinates (v, r). The function fi(v) 
encodes the fluctuations of dimensionless amplitude /io- The step function $(v ) in relation 
( |2.3|) indicates that the black hole results from the gravitational collapse of a massive (with 
mass M) null shell propagating along v = 0. Therefore inside the collapsing null shell the 
spacetime is flat. 

The conformal diagram of the whole geometry ( |2.1| ) in the absence of fluctuations 
(that is for /i = 0) is schematically shown in Figure p]. The dashed line on this figure 
shows a radial null ray which reaches J + at the moment of the retarded time u and which 
was sent from J~ at v = V(u) of advanced time. In the presence of the fluctuations, the 
conformal structure of the spacetime remains the same, but the function V(u) is modified. 
We shall study this modification in Section 4. 

Upon substituting the metric ( j2.1|) into Einstein's equations, one easily gets that the 
right-hand side of these equations takes the form 

T,j,v = [M5(v) + M(jl{v) #(«)] I, l v , (2.5) 



3 



r=0 




Figure 1: Conformal diagram for a black hole created by a collapse of a massive null shell. 
Solid dark line v — represents the collapsing massive null shell. 



where = —v yfl is a future-directed null vector tangent to a radial in-going null ray. 
From a "classical" viewpoint, n(v)l^l u is a fluctuating flux of positive (for 27m < uv < 
27r(n+l)) and negative (for 27r(n + l) < uv < 27r(n+2)) energy density. Correspondingly, 
the position of the apparent horizon r AH = 2M(1 + fi(v)) also fluctuates near its average 
value f AH = 2M. 

Following York |21[] we assume that the dimensionless amplitude /i = a(m Planck /M), 
where a is a pure number. In particular this assumption means that fiQ <C 1 for black 
holes of mass M ^> m Planck . It is also is assumed that in order to get a more realistic 
result one should average over a spectrum of metric fluctuations. 



3 Radial Ray Propagation 

3.1 Radial null rays in perturbed geometry 

We first study the propagation of radial null rays in the fluctuating black hole geometry 
Q2.1|) since we shall use, as usually done, geometric optics to construct the solutions of 
the wave equation. 

In-going rays are given by v = const and out-going rays obey the equation 

Adv = 2dr. (3.1) 
In order to solve this equation, we use a method of perturbations and write 

r = r(v) = R(v) + p(v) + a(v) + . . . = 2M \r(v) + p{v) + a(v)} + ... . (3.2) 
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R(v) is the solution of equation ( |3.1| ) in the absence of fluctuations, and p{y) and a{v) 
are respectively the first and second order perturbation in /i . Higher order corrections 
are denoted by dots. In what follows, we shall also often use the dimensionless versions 
of R, p and a which we mark by a tilde. 

The equation for out-going rays in the unperturbed metric, [ ( )' = d/dv ] 

A- |(i-3^4) 

can be easily integrated. Let us choose the value of the retarded time u and denote by 
r = R(v; u) the unperturbed trajectory of a radial ray which arrives to J + at the chosen 
time u. This trajectory can be found by solving the equation 

u = v — 2R* = const . (3.4) 

Here 

Rt = R ~ 2M + 2M In - (3.5) 

is the usual tortoise radial coordinate. 

The equations for the perturbations p(v) and a(v) are obtained by linearizing ( |3.1| ). 
Both functions obey the same equation 

M 

for v > and 

/ = 0, (3.7) 
for v < 0. For the first order perturbation, one has 

M , , 

f = P, F = -^' ( 3 - 8 ) 

and for the second order perturbation 

„ M M , . . 

f = °> F = lpW-R? p - (3 ' 9) 

In these equations, the retarded time u is a fixed parameter which specifies the unper- 
turbed ray under consideration and R = R(v; u). 



3.2 Perturbed horizon 

Before giving the general solution of the equations for the perturbations p(v) and a(v) 
we discuss the particular solution which describes the event horizon in the fluctuating 
geometry. 

First notice that R = 2M satisfies the unperturbed equation ( |3.3| ). This degenerate 
solution describes an outgoing null ray propagating along the unperturbed event horizon. 
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Starting with this solution we easily obtain the following solutions for the dimensionless 
perturbations 

Q cos(Qv) + sin(fH>) 



/'0 



Peh '■ 
2fl 2 (2 



l + tt 2 



fl 2 ) cos(2fi£) + Q(l - 5Q 2 ) sm(2Qd) 



2(l + Q 2 ) 2 (l + 4fi 2 ) 



where 



KV . 



(3.10) 
(3.11) 

(3.12) 



and k = (4M) _1 is the unperturbed surface gravity of the black hole. 

For both solutions it is possible to add a solution of the homogeneous equation. Such a 
solution corresponds to radial null rays propagating near the horizon in the unperturbed 
geometry. For these solutions the absolute value of p and a is infinitely growing. We 
choose the integration constants to exclude these solutions so that the perturbed radial 
null ray neither goes to infinity nor to r = 0. Therefore it describes the position of 
the event horizon in the perturbed geometry. It is easy to verify that ( |3.10|) and (|3.11 



coincides with the solution obtained by York [21] [see equations (4.8) and (4.9) of York's 
paper] . 

It is also interesting to compute the modified value of the surface area A of the event 
horizon. When averaging over time v, we find 



A = 47r(r h 2 or ( 



v)) = 16vrM z 



1 + 



/'o 



2(l + ft 2 )_ 

Similarly the average value of the surface gravity in the fluctuating geometry is 



(3.13) 



K 




1 + 3(p 2 ) - 2/i (psin^i;)) 



1 + 



2(1 + n 2 ) 



(3.14) 



Upon computing the modifications of the Hawking flux, we shall see that this 'renormal- 
ized' surface gravity will determined the modified temperature Tjj = R/2tt. The change 
of the area, 5 A = A — A, and the change of temperature, 5Th = Th — Th, of the black 
hole induced by the metric fluctuations obey the relationP] 

5 A _ 6T H 
~A ~ ~T^~ ' 



(3.15) 



3.3 Perturbed radial rays 

We now consider the general case, that is we assume that R ^ 2M. Equation 
can be easily integrated and gives / = const everywhere inside the collapsing shell. To 



It is interesting to note that similar corrections to the black hole surface area and temperature were 
obtained by Hu and Shiokawa p3| in their stationary model of metric stochasticity. The relation between 
5 A and STh in their case contains an additional factor 2. 
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integrate equation (|3.6| ), we change variable v to R(v;u), where as earlier the retarded 
time parameter u is fixed. This allows us to rewrite (|3.6|) as 



2M 
~R 



dR R? J 



(3.16) 



The solution of this equation is 



/ 



2M 



IF 



R 



} 2M\ 



■dR+f 



(3.17) 



The integration constant fo corresponds to the possibility to add a solution of the homoge- 
neous equation; that is, a solution of the equation (|3.16p with F = 0. We put f = since 
this choice is equivalent to the requirement that the ray propagating in the "perturbed" 
geometry arrives to J + at the same time u as the "unperturbed" ray r = R(v; u). 
In what follows it is convenient to introduce the dimensionless quantities 



R-2M 



x 



x 



2M 



1 1 



KU 



R* 



x + In x , f 



u 
AM 

f 



(3.18) 



(3.19) 



2M ' J 2M 

In these notations solutions (|3 . 1 7|) for the dimensionless perturbations p and a take the 
form 



p(x) 



X 



1 + x 



-Kx) 



!{*) = f§ (1 + MO 



X 



P(0 



The fluctuating mass term /i(^) = p(v(^)) which enters these equations is 

/i(0 = A* sin[n(f + lnf + 



(3.20) 
(3.21) 

(3.22) 
(3.23) 



To study quantum black hole radiation in the fluctuating geometry, we need to solve 
the wave equation in this geometry. As usual, we use the geometrical optic approxima- 
tionF]. Thus we only need to solve the following problem ]T|]. Consider a radial ray, which 
leaves 3~ at some advance time v and reaches J + at some retarded time u (see Fig. 1). 
What is the relation between u and v ? To establish the relation v = V(u) we use the 
above solution for a ray propagating in the fluctuating geometry outside the collapsing 
massive null shell, and glue it to the solution inside the shell. The latter means that the 



2 The validity of this approximation follows from the fact that the initial frequencies involved in the 
processes occuring at large u times are much larger than the characteristic frequency k — 1/4M. 
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values of r coordinate for both rays must be the same on the shell v = 0. Using this 
condition and the reflection at r = 0, one finds for V = V/4M 

V = - [1 + x + p(x ) + a(x )\ , (3.24) 

where xq is the value of x on the massive null shell (i.e. at v — 0). It depends on u and 
is a solution of the following equation 

u 4~ Xq -\- In xq = . (3.25) 

We shall be interested in the null rays arriving to J + at late time u. For such rays the 
equation (|3.25| ) can be solved by iterations and gives 

Xo = e" s (1 - e _fi ) + 0(e" 3S ) . (3.26) 

This relation shows that for the late-time regime Xq is very small. 



4 Calculation of V(?x)-function 
4.1 First order corrections 

In this section we analyse the perturbations p(x) and &(x) in the late-time regime. Before 
proceeding to the computation of these perturbations it is appropriate to make a few 
remarks. First, we notice that it is necessary to compute a, the quadratic fluctuation in 
/io, since our aim is to obtain all quadratic corrections to the flux and to the asymptotic 
spectrum. Secondly, it should be noticed that the following developments for determining 
V (u) in the perturbed geometry are quite similar to those of refs. |24], |25| which concerned 
the determination of V(u) when the energy-momentum tensor of Hawking quanta is taken 
into account. 

We start with the first order perturbation p(x), see eq. ( |3. 2Q|) . To calculate I(x) we 
notice that ( [Pi] ) can be written in the form 



where for n > 1 



and 



I(x) = po Im [PU^x- Q) + PU 2 {x- Q)} , (4.1) 
PU n (x;Q) = e inu P n (x-n) (4.2) 



e 

e 

By integrating by parts, it is easy to show that 



P n (x-Q)= I (4.3) 



[n — 1 — ill) x n 1 IU n — 1 — ill 
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Using this relation, we can rewrite (|4.1|) as 

I(x) = fi Im {-^— \ e «K*+*+**)l + PUl {x- Q)} } 
L 1 — lit L x J J 

sin[f2(£i + x + In a;) + arctanfi] + /(x) , 



/'o 



where 



/(x) = /i Im 



(4.5) 



(4.6) 



i -in 

Function Pi(x;Q) can be expressed in terms of the incomplete gamma function T(a,()- 
Using the relation 

l ^d^ = fi- a T(a,fiu), (4.7) 



we get 



Px{x; G) = {-iO)- iU r(zO, -iQx) = e -^-^^ r (ifi, -ifix) . (4.8) 
Furthermore, the incomplete gamma function T(a, () allows the following series expansion 

C" 



T(a,0 = r(a) 
Thus for small x (fix <C 1) we have 



n=0 



T(a + n + 1) 



(4.9) 



T{ifl, -iflx) « r(ifi) + ^ e 2 ? x iQ e ininn 



and 



Pi(x;0) « e 



r(ifi) + 



(4.10) 



(4.11) 



For the calculation of V(w) we need to know I(xq) where x is the solution of equation 
( ^251) . Using (|J), one gets 



I(x ) 



x (i + n 2 ) 



+ iM , 



where 



I(x Q ) = no Im 



te 



+ /iQ Im 



e 2 



-ifi In f2 



-mo r(tfi) 



(4.12) 



(4.13) 



For small xo the first term in the right-hand side gives [io/[Q(l + fi 2 )]. To calculate the 
second term contribution, we notice that 



e 2 



where 



-ifl In f2 



r(zO) = q(tt) e 



-j(/3 r (0) 



2^ 



(4.14) 



(4.15) 
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and <pr{ty is a real function which for large values of Q, Q — > oo, has the following 
expansion 



1 



-lY^B 



2/i 



^ ) "° + 4 7r + £(2n-l)(2n)^- 
where B n are Bernoulli numbers. 

Inserting ( 4.14 ) into (|4.13|) one gets for small Xq 



(4.16) 



Uxo) 



Ik) 



0(1 + n 2 ) ru VTTf 



sin $i 



where 



$i = $(fi) + arctanfi, $(0) = fiw - y? r (^) • 
Collecting all the results, we finally obtain 

n i q(n) 



I(X ) « yU 



a;o(l + fi 2 ) fi(l + fi 2 ) x/TTTF 



sin $i 



(4.17) 
(4.18) 

(4.19) 



This result can be used to obtain the asymptotic form of the function V(u). Using 
equations (|3.20|) , (|3.24|) , (|4.19|) , and (|4.20|) , we find up to the first order in e~ u 



- V « V + e-* [C x + C 2 sin($i)] + V 2 (u) 



where 



V = 1 + no 



n 



i + n 2 



6-1 a C - a ^ 

^l — ^—^o 777- — ; 09 n j ^2-/^0 



and 



n(i + n 2 ) 

V 2 (u) = a(x (u)) 



(4.20) 

(4.21) 
(4.22) 



is the second order correction. 



4.2 Second order corrections 



To calculate the second order correction a, we must estimate the integrals of eq. (|3.22|) 
By integrating by parts, a can be rewritten as 



x 



I 2 (x) 



aix) 



1 + x \ 2(1 + x) 2 



where as earlier 



and 



It is easy to see that 



I(x) = h(x) + I 2 (x) , 



(4.23) 

(4.24) 
(4.25) 

(4.26) 
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Thus we have 



x 



a(x. 



1 + x 



where 



and 



T(x) =/i Im [e iQil Q 2 (x;n) 



(4.27) 

(4.28) 
(4.29) 



The functions I n (x) are related to PU n (x; Q) defined by equation Q4.2Q as follows 

I n (x)=iM)hn.[PU n (x',n)]. (4.30) 

Using the results of the previous section, we can obtain the expansions of I n (xo) for small 
Xq. Thus, the problem of finding a(xo) for small xo is reduced to study the function T(x). 
For this purpose we first obtain a recursion relation for Q n (x; O). By integrating by parts, 
we have 



Q n (x\ Q) 



e inx h(x) 



1 



n 



1 - in) x™- 1 -^ (n - 1 - ifl) A f ' 

in 



n— if2 



+ 



■Q n _i(a;;fi) 



, (4.31) 

n — \ — i\l 

Using ( |3.23| ) and the definition ( |4.2j ) of PU n (x; fl) we can write the second term of (|4.31|) 

as 

rno AC J®£ T,(C\ 

PU n (x- 2ft) 



— ^In — - _ 2^° 



1 



(n — 1) x 



n— 1 



(4.32) 



Relations (|4.31| ) and (|4.32|) allow us to write Q 2 {x] ft) as 



Q 2 (x; n) 



/i e 



(l-ift)x 1 -^ 2i(l-in) 

in 



PU 2 {x;2Q) - 



x 



+ 



We also have 



Qi(x; n) = 



2i 



DO 



d£ e 



1 - in 



n) . 



(4.33) 



iQu 



€ ?7 



1-ifi 



d£ e iQ t f°° drj e 
1 



1-iQ 



-iQrj 



- ^2 p-ifiS 

" 2z 



- (P£Mx;ft)) 2 -S(x;ft) 



where 



(4.34) 
(4.35) 



This function allows the following representation (see Appendix A) 

5(ar; n) = ~ e- 2 ? 4 "^ 111 " r(-in) x iQ iF x (in; 1 + ifi; ifix) 
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+ - [\n(iQx) + f(-iQ; itix) - ^(ifi)] . 
where iF\ is a hypergeometric function, ip(z) = (dlnT(z)/dz), and 



f(a;y) = ^ 



y 



Here 



[1 + a) 



n(l + a) n 
r(l + a + n) 



r(i + «) 

Combining these results, we can rewrite (|4.33|) as 



(4.36) 



(4.37) 



(4.38) 



e lUu Q 2 (x; fi) 



1 — zO) x 



2i(1 -ift) 



P?7 2 (a:;2fi) 



+ 



(4.39) 



2(1 -ifl) 

The final expression for the second order perturbation a(xo) is obtained by inserting 
flQgp , flC59D , and flOOD into (|Q7D . 

Since we are interested in 5"(:co) where Xo is small we need only to know the first terms 
of its expansion in powers of Xq. To do this we use the following asymptotics at small 
value of Xq 

(4.40) 



PU l {x Q -n)^- + q{n)e^ 



PU 2 (x ;n) 



in 



i 

Xq 



iSlPU x {x Q \Sl) 



S(x ; n) « — 



171 



(4.41) 



(4.42) 



The approximative values of Ii(xq) and ^2(^0) can be easily obtained by using the relation 
( f4.30[ ). The calculations the asymptotic value of u{xq) is hence straightforward but quite 
long. We use Maple to perform these calculations. The result is 



a(x ) « $ [cr + (TiXq + a'} 



where 



^0 



Q 2 (2-Q 2 



(l + fi 2 ) 2 (l + 4ft 2 ) ' 
a 1 = a° 1 + a\ q(0) + a\ q 2 (Q) + a\ q(2Q) 
_ 2 + 3Q 2 - 33ft 4 + 20ft 6 

Gl ~ 4ft 2 (i + ft 2 ) 2 (i + 4ft 2 ) 

(ft 2 - l)[ftcos($(ft)) + sin($(ft))] 

" ft(i + n 2 ) 2 ' 

;i - ft 2 ) + [ftsin(2$(ft)) - cos(2$(ft))] 
4(1 + ft 2 ) 



(4.43) 

(4.44) 
(4.45) 
(4.46) 

(4.47) 

(4.48) 



12 



and 



a 



Q [(2n 2 - 1) sin($(2fi)) - 3ft cos($(2fi))] 
(l + ft 2 )(l+4ft 2 ) ' 

[vrft - 2 In(fixo) - 2Im[(fi - 



x 



(4.49) 



(4.50) 



4(1 + n 2 ) 

Let us recall that in these expressions Xq is a function of the retarded time u given 
by ( |3.26| ). We can now include the second order corrections (|4.22 ) into the expression for 
V(u) and write it in a form similar to ( J4.20p 



V « V + e- s [C a + C 2 sin($i) + C 3 sin($ 2 ) + C 4 sin($ 3 ) + tfu] 



where 



V = 1 + no 



n 



l + tt 2 



Q 2 (2-Q 2 ) 



d = l- no 



Q 2 -l 
Q{l + Q 2 ) 



1 + ft 2 ) 2 (l + 4ft 2 ) ' 

2 + 3ft 2 - 33ft 4 + 20ft 6 
4ft 2 (l + ft 2 ) 2 (l + 4ft 2 ) 



4(1 + ft 2 ) 



ft 2 ^ 2 (7rft-2mft-2Im[(ft-i)V>(ift)]) 



C*2 — /^O 



ff(n) 



x/i + ft 2 



4(i + ft 2 ) 
ft 2 -i 



ft(l + ft 2 ) 



°3 — A'o 



2 



C 



4 v a + n 2 

ftg(2ft) 



l + ft 2 )(l +4ft 2 ) 



2(1 + ft 2 ) 



and the phases $j are defined as follows 



$! = Qu — y?r(ft) + arctanft , 

$ 2 = 2Qu - 2y? r (ft) - arctan(^) 

3ft 

$ 3 = 2Vtu - y> r (2ft) - arctan( 



(4.51) 
(4.52) 



(4.53) 
(4.54) 
(4.55) 
(4.56) 

(4.57) 

(4.58) 
(4.59) 
(4.60) 

'2ft 2 -l y v ; 

Notice, that the coefficients Vo, C\ and Ci which appeared in the first order expression 
( f4.20| ) are now replaced in ( |4.51| ) by the new coefficients Vo, C\ and C2 with the only 
difference that the corresponding coefficients get second order corrections. Notice also 
that in (|4.51|) , the terms with double frequency 2ft and a term which is linear in u are 
new with respect to the first oder result (|4.20|) . 
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5 Calculation of the Flux of Hawking Radiation 



Now we derive the s-mode contribution to Hawking radiation. In what follows, we shall 
neglect the scattering by the gravitational potential barrier which appears in the 4D 
Dalembertian. In other words, we use 2D approximation in which ingoing and outgoing 
modes completely decouple. This strong hypothesis requires some explanations. The 
decoupling of the modes greatly simplifies the calculation of the asymptotic flux when the 
metric is no longer static. Indeed, the height of the potential barrier now depends on v 
in the metric eq. ( j2.1|) . Therefore one looses the fact that the transmission coefficients 
are diagonal in energy. Moreover, the new coefficients will also mix positive and negative 
frequency modes. This will lead to additional pair creation probabilities. Thus there 
will be interference effects between the usual pair creation amplitudes induced by the 
frequency mixing governed by eq. ( [4.51] ) and these new coefficients. 

To determine the importance of these interesting effects is complicated and goes be- 
yond the scope of the present paper which is to describe the effects on Hawking radiation 
induced by the fluctuations of the geometry in the very close vicinity of the horizon. In 
this respect, we wish to emphasize that our classical metric has been choosen to mimic 
the near horizon quantum fluctuations and not the fluctuations of the height of the bar- 
rier around r = 3M. On physical grounds, in a self-consistent treatment, one might 
expect that the residual fluctuations around 3M be much smaller that the near horizon 
ones. Therefore the neglection of the time dependence of the barrier might turn out to 
be physically legitimate. 

In the 2 dimensional simplified description, when the field is in its vacuum state before 
the formation of the black hole, the mean energy flux at J + is 



dE 
du 





-1/2- 



(5-1) 



Here V(u) is the function calculated in the previous section. Notice also, that u in this 
relation is the proper time at J + in the perturbed geometry, see the remark made after 
eq. ( p. 17] ). Thus this is the time which defines positive frequency at J + . 

Before presenting the results of the calculations of dE/du we make several remarks. 
First, it is evident that the expression for dE/du does not depend on the value of constant 
Vo in equation ( |4. 51|) . For this reason we can put it equal to zero. This corresponds to a 
simple redefinition V — > V + Vq. Moreover, dE/du is not changed if we multiply V by an 
arbitrary constant. For these reasons the calculations of dE/du can be performed with 
the simplified form for V 



V = -e 5 [1 + Ai sin(fiw + y?i) + A 2 sin(2fiw + y? 2 ) + Cu] . (5.2) 
In this expression, we have kept all terms up to second order in /xq and introduced the 
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following notations 

Ai = fi ai + filh 



Ho 



l + 2/xc 



Q 2 



1 



Q(l + tt 2 ) 



(5.3) 



ipx = — (fr + arctanfi . (5-4) 

In the same way the two terms in ( f4.51| ) with coefficients C3 and C4 having the same 
dependence 2Qu on the retarded time u have been combined into the following single 
term 

A 2 sin(2fiw + (p 2 ) = C 3 sin($ 2 ) + C 4 sin($ 3 ) (5.5) 

A 2 is of second order in fi and, with notations to ( |5.3|) , can be written as A 2 = nl\ b 2 . The 
explicit expressions for b 2 and ip 2 can be obtained easily, but since they are very long and 
are not important for our final result we do not reproduce them. Note also, that since C 
is already a quantity of second order in //q, the redefinition of V does not affect it. 

The calculation of dE/du from (|5.1| ) and (|5~2|) is straightforward but long and was 
performed by using Maple. It is convenient to write the result in the form 



dE/du = (dE/du) pcim + (dE/du) 



(5.6) 



where (dE/du) pcrm is the mean value of the flux and (dE/du) auct is its fluctuating part. 
The latter will of course not contribute to the total energy received on J + . 
The constant part is 

(dE/duY™ 



and the fluctuating part is 

(dE/du) 



48tt 



-yfiVl + n 2 q(n) cos(ftw + ip x ) 



1 



l + ^0 2 g 2 (fi) -2C 



(5.7) 



l-VL 2 VL 2 (1 
q(fl) = cos(fH + (fx) + q 2 (Q) 



5Q 2 



+ q 2 (n) 



Q(l + 4Q 2 ] 



8(1 + Q 2 ) 



cos(2fiw + 2(p\ 



sin(2fiw + 2(pt) - b 2 Vt(l + 4fi 2 ) cos(2fiw + Lp 2 ) 



(5.1 



4(1 + n 2 ) 

The function q(Q) is given by ( |4.15|) and, in the last term, b 2 cos(2fiw + <p 2 ) is equal to 



b 2 cos(2f>u + <p 2 ) 



q 2 (ty 

4 V / TTT^ 



cos($ 2 ) 



ttq(2tt) 



(1 + Q 2 )(1 + AQ 2 ) 



■ cos($ 3 ) 



(5.9) 



The remarkable fact is that, to second order in fi Q , the correction term which is linear 
in u in Q5.2D does not give any time-dependent contribution. It only gives an additional 
constant to (dE/du) pcTU1 in ( |5.7|) . This has the following simple explanation: the term Cu 
can be removed from ( |5.2|) by absorbing it into e~ u without changing the other terms in 
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our second order expression. This transformation corresponds to the "renormalization" 
of the surface gravity 



k{1-C) = K 



1 + 



/'o 



2(1 + fi 2 )_ 

Hence the expression for (dE/du) pciin can be identically rewritten as 



(dE/du) pc 



-2 r 



48tt 



n--rtjnV(n) 



48vr 



1 + 



exp(27rf2) 



(5.10) 



(5.11) 



It is interesting to note that the renormalized surface gravity k t which is introduced here 
coincides with the average value of the surface gravity R defined by equation ( |3.14| ). 



6 The Modified Asymptotic Spectrum 

Instead of focusing on integrated quantities, we shall now consider how the asymptotic 
spectrum is modified by the fluctuating part of the metric. As usual, the asymptotic 
spectrum is characterized by the angular momentum (that we take to zero) and the 
asymptotic energy, A, the eigenvalue of id u . Indeed the fluctuations we are considering 
do not affect the stationary character of the asymptotic (large r) metric. 

To obtain the modified spectrum, we need to compute the Bogoliubov coefficients in 
the modified geometry. Before proceeding to this calculation, it should be noticed that the 
2D expressions we shall use are exact only for large A, i.e. n\ ^> 1. At lower frequencies, 
there is indeed a potential barrier in the 4D Dalembertian for s-waves which reduces the 
transmitted flux in a static geometry, c.f. the discussion at the beginning of the former 
Section. 

We first recall how the well known properties of Hawking spectrum are extracted from 
the Bogoliubov coefficients. The latter are given by the overlap of the initial (infalling) 
modes which are specified on J~ and the final (outgoing) modes specified on J~ + . Both 
are solutions of the Dalembertian equation in the metric ( |2.1| ). For s-waves and under the 
neglection of the potential barrier, these modes satisfy the 2D equation d u d v <p = 0. Thus 
the in-modes can be decomposed in terms of planes waves 



<t>v(v) 



V47TZ/ 

where v is the energy measured on J~ . Similarly, the out-modes are 



(6.1) 



-iXu 



0a(«) 



(6.2) 



where A is the energy measured on J + . 

The scattering of in-modes in the time dependent geometry simply follows from the 
'reflection' condition on r = wherein the Wronskian must vanish. This implies that the 
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scattered in- modes are given by 4> u (V(u)). Then the Bogoliubov coefficients are 

p <_» r e ivV(u) e ~i\u 

a u ,\ = du (j)* u (V(u)) id u (p x (u) = du 



i 



p <_» p e ~iuV(u) e ~i\u 

Pu,\ = du 4> u (V(u)) id u 4> x {u) = du — 7 == — r — _ (6.3) 

The second expressions follow from integration by part and the neglect of the end-point 
contributions. 



In the unperturbed geometry, for large u, kV{u) is given by —1 — e~ KU , see eq. ( |3.24| ). 
The validity of this asymptotic expression requires that the initial frequencies v be large 
enough. Indeed, upon applying the stationary phase condition to the integrand of the 
coefficient a u> \ one obtains the following relation for the position of the saddle point u sp 

\ = ve -^v . ( 6 .4) 

This equation simply specifies the value of u at which the Doppler shifted initial frequency 
v resonates with A. Notice in passing that it is this exponentially small Doppler factor 
which leads to the necessity of considering arbitrary large initial ('transplanckian') fre- 
quencies. From eq. (|6.4D, one deduces that large v means that the corresponding value of 
e -K,u sp sa ti s fi es e -™ S p ^ e ~2Ku sp ^ N ti ce a l S o that the stationary phase condition applied 
to the (3 coefficient, leads to the same condition up to an overall negative sign. This 
implies that the location of the saddle point in the complex u-plane receives an imagi- 
nary contribution equal to —itt/k. The determination of the sign of this imaginary part 
follows from the fact that V(u) appearing in e~ wV must belong to the lower half complex 
\^-plane. Physically this amounts to specify that the in-vacuum contains no excitation 
characterized by positive v. Mathematically it gives I/^a/c^aI 2 = e~ 2wX / K . 

These considerations based on a saddle point analysis are confirmed by the exact 
integration of eqs. (|6.3D. It will be found useful to express the exact expressions in terms 
of the following function 



B(y,X) = r(<-) Jt^-A-)- iX/k (6.5) 

where T(x) is the Euler complete gamma function. The norm of this function gives: 
\B\ 2 = (e 27rA / re — l) -1 . Upon extending the domain of validity of the asymptotic behaviour 
of V(u) for all u, one finds 




A,, A = B{v,\) x \— -e 1 "" 1 . (6.6) 



Two crucial properties follow from eq. (|6.6|) : firstly, the Planck distribution charac- 
terizing the mean number of out-quanta in the in-vacuum (up to a normalisation factor, 
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it is obtained from |/3^a| 2 ), and secondly, the existence of a constant flux of out-quanta. 
The stationarity follows from the fact that the phases of (3 and a are both proportional 
to y- %x l K , This implies indeed that the value of the energy flux is constant. To prove it, 
we recall that the renormalized value of the energy flux in the in-vacuum, when expressed 
in terms of the Bogoliubov coefficients, is 



dE 



r°° v'AA' 
d \ / d\'^± 

du Jo Jo 2ir 



-i{\-\')u 



(poo 



(6.7) 



By using eq. (|6.6| ) and performing the integral over u, one immediately obtains that the 
second term, the interfering one, vanishes and that the first term is constant. 

We shall now determine how these properties are affected by the fluctuations of the 



geometry described by eq. fl2.1|) . Under our restriction to s- waves and neglection of the 
potential barrier, it suffices to repeat the same procedure with the modified function V(u), 
as if we were considering 2D propagation. The asymptotic expression we shall use is given 
by 

- K r V(u) =V + e~ KrU [l + A x sin(cuM + (p x ) + A 2 sin(2u;w + tp 2 )\ , (6.8) 

where the values of V , A%, and A 2 can be found in eqs. ( |4.52| ), ( |5.3| ), and ( |5.5| ). We have 
kept the constant shift in V since it will modify the absolute phase of the Bogoliubov 
coefficients. The symbol n r designates the 'renormalized' surface gravity introduced in 
eq. ( |5.10| ). We have indeed absorbed the linear term in u of eq. ( |5.2| ) in the redefinition 
of k. As before, since n r — k = 0(/ig) and since we shall work up to quadratic corrections 
in /io, the modification of the terms proportional to A\ and A 2 are irrelevant. 

To reveal the nature of the modifications induced by the time dependent fluctuations 
of the metric, it is appropriate to analyse the (3 coefficient. Up to quadratic order in the 
metric fluctuations of amplitude /io, the modified (3 coefficient reads 



nmod 
Pu,X 



OO g 

du—j= 

-oo V7rA _1 
e -i\u e 



iXu p —ivV(u) 



du 



V47TZ/ 



-1 



2 Kj T 



l-rf— 



) 2 e- 2 ^ u A 2 1 sm 2 {uju + ! f l ) 



v 



+i—e KrU {A x sin(u;M + + A 2 sin(2cjw + ip 2 )} 



(6.9) 



By decomposing the sines into imaginary exponentials and by using several times the 
relation T(x + 1) = xT(x) we get 



4 



mod 
v.X 



1 



JuV /K r 



x \B r {v,X) 



A 2 

1 - — ^-(A - in r )X 



2wn r v 

\B r (v, A - lu)V\ -u e iipl - B r (u, A + + u e"^ 1 

2 Kj r 
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-V\— \B r (u, X - 2lu)V\ - 2uo e^ 2 - B r (v, A + 2u)^/\ + 2u e 

+V\~^ BJu, A - 2u) (\-2w- iK r WX-2u e 2l(pi 

+B r (y, A + 2u){\ + 2uj- iK r )VX + 2u e ' 2i<pi ] } . 



-up 2 



(6.10) 



Here B r (v,X) designates the function B(u, A) defined in eq. (6.5) with the renormalized 
surface gravity. 

The replacement of B(u, A) by B{y, A ± in the last three terms of ( 6.10Q indicates 
that u, the frequency of the fluctuating metric, enters into the expressions in such a way 
that the 'effective' frequency, i.e. the one which weighs the new amplitudes, is A ± uj. 

Physically, this leads to a modification of the mean (quantum averaged and time 
averaged) number of quanta which reach J + per unit u time. This averaged flux is 



mod 1 2 



(*A> 



(6.11) 



We have implemented time average by integrating over v up to the cut-off frequency iV 
and dividing the resulting expression by / dv j 1 n r v = Au. This last equality follows from 
the resonance condition, eq. Q6.4Q , and its validity requires that Au ^> l/«: r , as in usual 
Golden Rule estimates. 

In this ratio, to quadratic order in /xq, only four terms contribute. The first two terms 
of the first line of eq. fl6.10|) and those resulting from the square of the second line. This 
is because the integral of all other terms are crossed terms of the form B(v, X)B*(u, X±uj) 
which lead to oscillatory integrands containing e ± % ^ u sv( v ) _ These terms do not contribute 
to the production rate, eq. ( |6.11|) , in the limit of large kAu. 

By combining the four non-vanishing contributions we get 



(nx) 



2n [ e 2nX / Kr - 1 

2 hh r 



1 



X — UJ 



+ 



X + u) 



2tx(\—uj) / K r ^ g27r(A+o;)/K r J 



(6.12) 



Having obtained the spectrum of the particle number flux, the spectrum of the energy 
flux is 

dE , , . k t 



where A 



and 



du dX 

X/ K r , = Oj/K r , 



X {fix) 



2ty 



/(A) 



/(A) + ^F(A;fi) 



A 



exp(27rA) 



F{A;fl) 



ttA 2 



/(o) [/(a - n) + /(A + n) - 2/(A)] 



n 2 (i + n 2 ) 

Plots of the function F(A; Q) for different values of Q are given in Figure 2. 



(6.13) 

(6.14) 
(6.15) 
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Figure 2: Function F(A, Q) for different values of dimensionless frequency Q: Q = 0.01 - 
solid line, Q = 0.1 - dashed line, and Q = 0.5 - dotted line. 



In order to obtain the constant part of the density of energy flux one must integrate 
( |6.13|) over the frequency A 



dE 
du 



2^ 



dA 



f(A) + ^F(A;Q) 



(6.16) 



Using formulas ( P-l|) and ( |B.10| ) from appendix B, it is easy to verify that ( |6.16| ) coincides 
with the expression ( |5.11| ) for (dE/du) pcim we obtained earlier. 

The other physical consequence of the corrections terms to a and (3 is the following. 
The mean instantaneous energy flux now oscillates around its time average with the 
harmonics of to. This can be seens from the oscillatory terms B(u, X)B*(u, X ± (2)u) 
which behave in e ^( 2 )^ u s P M w h en parametrized by nu sp = lnz//A, the location of the 
dominant contribution to a v ^\. The amplitudes of the fluctuations are linear in A\ and 
A 2 . Moreover, the phase shifts of the modified coefficients with respect the unperturbed 
ones differ for a and (3. This implies that the second term in eq. ( |6.7| ) no longer vanishes 
even though it still does not contribute to the mean (time average) energy flux. 
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7 Conclusion 



In this paper we studied how the fluctuations of the geometry near the black-hole horizon 
affect Hawking radiation. To characterize the fluctuations of geometry, we used York's 



model [[H]] in which they are described by an in-coming Vaydia metric with a time depen- 
dent mass m(v) which fluctuates around its mean with amplitude fio and period 2ir/u. 
For further simplicity, we only considered the s-modes of a quantum scalar massless field 
and we also neglected the scattering by the gravitational potential barrier. Using these 
simplifying assumptions, we reached the following conclusions. 

First, the expectation value of the outgoing flux of energy is no longer constant. It 
now fluctuates around its time averaged value with frequencies given by harmonics of u 
and with amplitudes starting with a term linear in /i , see eq. ( |5.8| ). The fact that the 
phase of the fluctuations of the expectation value of dE/du is well defined results from 
the fact that the fluctuations we considered were treated classically. In a more quantum 
mechanical treatment, these well defined phases will probably be replaced by a more 
diffuse ensemble of phases. 

Secondly, the time averaged value of the outgoing flux of energy is modified. One 
part of this modification is connected with the renormalization of the surface gravity of 
the fluctuating black hole given by expression ( |5.10| ). The other part is an additional 
factor given by yUQ7rf2/(exp(27rf2) — 1), see eq. ( |5.11| ). Both changes are second order in 
yUo- More surprisingly, they decrease for large Q. Indeed, one might have feared that fast 
metric fluctuations would lead to copious pair creation 0. We conjecture that we do not 
find abundant production since we worked in a 2 dimensional model in which the time 
dependence of the metric does not affect directly massless modes, thanks to conformal 
invariance. 

Thirdly, the asymptotic spectrum of Hawking radiation is also modified. Besides the 
renormalization of the surface gravity which shifts the temperature, the modified spectrum 
( |6.12|) contains three additional correction terms. The two last terms in that equation 
contain Bose thermal factors of the form 1 / (exp(27r(A ±u) /k) — 1). In these relations, the 
frequency of geometry fluctuations, ±u, plays the role of a chemical potential. The pres- 
ence of such chemical potential is reminiscent to superradiance. This fact supports the 
general ideas proposed by York since the appearance of these factors might be expected 
from the existence of a quantum ergoshpere. Indeed, due to quantum fluctuations, the 
average position of the event horizon is moved by a term proportional to the second power 
/!q of the amplitude of fluctuations, while the temporal position of the apparent horizon is 
fluctuating with amplitude /x - An alternative way to describe these fluctuations is to say 
that there exists a blurring of the physical null cone at the unperturbed horizon. Because 
of the existence of negative energy states inside the unperturbed black hole matter can 
escape from the narrow region close to the horizon. This leakage of energy is seen as 



3 We are grateful to T. Jacobson for pointing out this fact to us. 
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Hawking radiation ||21|| . Under the same conditions one can expect an additional ampli- 
fication of Hawking quanta while they are propagating close to the fluctuating horizon. 
The amplification factor we got in the expression for the modified spectrum of Hawking 
radiation may be considered as an indication to this effect. 

The modifications in the black hole temperature and surface area in the presence of 
metric fluctuations raise the question about the modifications of black hole thermodynam- 
ics. If we identify the energy of the system E with the averaged mass of the black hole 
M, and the temperature of the black hole, T, with K r /27r, then the first law dE = T dS 
and eqs.( P-13| ) and ( |5.10| ) define the averaged entropy to be 



Therefore, one looses the relationship between the entropy and a fourth of the area. 
If one writes the amplitude of the fluctuations as /io = am P i anck /M where a is dimen- 
sionless, one has 

M 2 

S = Att— s, (7.18) 

^^Planck 

where 

Area 2 , 
.= 1+^ (7-19) 

does not depend on the black hole mass. It is worth noticing that this modification of the 
black hole entropy is exactly of the same form as in theories with corrections quadratic 
in the curvature Wh - We recall that these corrections arise in the effective action by 



taking, in the one-loop approximation, the average of the gravitational equations over 
quantum fluctuations of the metric. This observation might give a possible explanation 
to the origin of the similarity between ( [7.18 ) and the results of Ref. [27]. Similarily, it 



would be interesting to find the relation between the results obtained in this paper and 
the quantum treatments of p9T[-||31 . 



Even though these results were obtained in an extremely simplified model in which 
the metric fluctuations were treated classically, we believe that they indicate what might 
be the impact of the quantum fluctuations of the near horizon geometry on black hole 
radiance. 
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A Calculation of Function S(x). 



To calculate function 



S(x) 



°° d£ ^ j°° drj e~ inri 



£1— if2 yyl+iQ 



= (in) in I 2^r(-in,ine), (A.i) 

Jx t; 

we use the following general result, which can be found in |26| (volume II relation 2.10.3.6) 

x™" 1 (x-af- l e cx T{u,cx) dx 

a a+f3-l T ( u j 1 _ a _^ iFi ( a] a + /3 . ac j 

B(P, 1 — ct — (3 — v) 2 F 2 (a + 1; v + 1, a + (3 + v\ ac) 



TCC 



l ~ a -P T(a + (3-l) 



Here 



sin[(a + p + u)ir] T(l - v) 
B(x,y) - 



iFi(l-/3;2-a;-/3; ac) . 



(A.2) 



(A.3) 



T(x + y) 

This relation is valid for a > 0, Re/3 > 0, Re (a + + v) < 2; | argc | < it. For a particular 
case (3 = 1 one has 

a a T(z/) 



x a " 1 e ca; r(z/,ca;) 



a 



1 + a; ac) 



+ 



i/(a + z/ 



2 F 2 (a + v, 1; 1 + z/, 1 + a + v\ ac) 



+ — 



TCC 



r(a) 



sin[(a + i/)7r] T(l — z/) 
The integral in (|A.l ) is of the form (|A.4 ) with 



a = iQ , z/ = — ifl , a = x, c = ifl . 



(A.4) 



(A.5) 



For this case a + v = 0, and the last two terms in ( |A.4|) are divergent. To consider 
this limit, we define 7 = a + v and rewrite the sum of these two divergent terms in the 
following form 



a a+v c v 
a(a + v) 



2 F 2 (a + v, 1; 1 + v, 1 + a + v\ ac) + 



TIC 



T(a) 



sin[(a + u)tc] T(l — v) 



^—7 



where 



z = I 

7 



^,F,( 7 ,l;l + „,l +r ,») + .V^lrl 

sm(77r) 1 (1 — z/) 



(A.6) 
(A.7) 



23 



and y = ac. 

Using the expression for the hypergeometric function 

oo n 

( 7 + n )(l + i/) T 

where (1 + z/) n = T(l + n + z/) /T(l + i/), we can write for small 7 

2 F 2 (7, 1; 1 + 1/, 1 + 7; y) « 1 + 7/(1/; y) , 



where 



We also have 



00 r. 

/(*;*) = £ 



^ n(l + l>) n ' 

r (7"^) ^ 1 , r'(-z/) 1 1 

7= n-fjT- 



r(i -1/) 1/ r(i -1/) ' 1/ 1/ 

where V^- 2 ) = <ilnr(z)/(iz. 

Substituting expression ( |A.9| ) and ( |A.11| ) into (|A.7|) and expanding sin(77r) 



| 7 2 tt 2 + ...), we get 

Z = -[^(1 + 7/07 2/)) -l-T^-f)] 

7 

- 1 



7 

Thus for 1/ = —a we have 



r°° a a T(— a) 

/ x a ~ l e cx T(-a, cx) dx = iF 1 (a; 1 + a; ac 

J a a 



a 



[ln(ac) + /(—a; ac) — ip(a)] 



Using this result, we obtain for S(x) the following expression 

S(x) = ^ e ininn+ininx e'^' 2 r(-ifl) iF^ifi; 1 + ifi; zftz) 

+ 1 [ln(iOa;) + /(-ifl; ifix) - ^(ifi)] . 

B Calculation of Integrals. 

In this appendix we demonstrate that 



/o exp[27r(A + ft)] -1 Vo exp[2-7r(A - Q)] — 1 120 
First we notice that 

J= lim [J + (p) + J_(p)] 
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where 



J±(P) 



°° dX A 2 (A±ft)exp(-p|A|) 



exp[2vr(A±ft)] - 1 

Making change of variable of integration A — ► —A in the expression for J-(p) we get 

o dX X 2 {X + Q) exp{-p\X\) 



(B-3) 



J_(p) 
Thus we have 
J+(p) + J_(p) 



dAA 2 (A + ft)expG»A) + 



dX A 2 (A + n) exp(pA) + 



3 exp[2vr(A + fi)] - 1 

00 dx (x — Q) 2 xexp(— p\x — Q|) 
-oo exp(27rx) — 1 



(B.4) 



(B.5) 



In the last integral in the right hand side we made change of variables A = x — Q. 
Consider now the integral 

r° dx (x — Q) 2 x exp(—p\x — Q\) 
J-oo exp(2-7rx) — 1 

By changing variables x — ► —x it can be identically rewritten as 

r°° dxx (x + Q) 2 exp(-px) 



(B.6) 



exp(— pQ) / dxx(x + Q) exp(— px) + exp(— pQ) 



exp(27nr) 



(B.7) 



Using ( B.5|) and ( B.71 ) and taking the limit p — > we get 

ro ° dxx (x 2 + fl 2 



J = 2 



o exp(27rx) 



+ AJ. 



(B.8) 



where 



/0 roo 
dxx 2 (x + Q) exp(px) + exp(— pQ) / dxx (x + ft) 2 exp(— px) 
-oo ' Jo 



lim 



' 6 2ft* 

p A p 3 



+ e 



-pQ 



' 6 4ft ft 2% 

^4 ^3 ^2 



-Q 4 



The integrals which enter ( |B.8D can be easily calculated 



dx x 



k exp(2vrx) - 1 24 
Combining these results we get (|B.1|). 



o exp(27nc) 



12 



1 

240 



(B.9) 



(B.10) 
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